Path Entropy Changes in Adiabatic Approximation 
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By applying adiabatic theorem to a Markovian system, we calculate the adiabatic and diabatic 
entropy changes along a path. As well known, the total path entropy change is separated into two 
parts, system and environment entropy changes, AStot = AS sy s + ASenv. The environment entropy 
change, ASenv, is divided again into two parts, an adiabatic contribution due to work, ASw, and 
a diabatic contributions due to heat, ASg. In an adiabatic process, total path entropy change is 
same with the adiabatic path entropy change, AS a, which is given by sum of system entropy change 
and adiabatic contribution, AS a = ASsys + ASw- Mathematical form of AS a is a type of excess 
heat entropy change, but ASa is due to work. By which, it is shown that the terms adiabatic and 
non-adiabatic contributions of AS na and ASa in [Phys. Rev. Lett. 104, 090601 (2010)] should be 
completely switched, i.e. AS na — > ASa and AS a — > ASq in fact. 

PACS numbers: 05.70.Ln, 05.40.-a, 02.50.-r 
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Introduction.— One of the attractions of fluctuation 
theorem (FT) is that they successfully describe the be- 
havior of the second law of thermodynamics (AS) > in 
microscopic limit. Ever since the first FT is introduced 
by Evans, Cohen, and Morriss [l| in 1993, various FTs 
have been developed. Among them the FTs proposed by 
Esposito and Van den Broeck [2j, yj are considered to be 
the most generalized ones. 

In their study, Esposito and Van den Broeck have 
expressed the total entropy change for a single path 
as a sum of two distinctive contributions, A.S na and 
AS a , which are called respectively as "non-adiabatic" 
and "adiabatic" because of the (quantum) adiabatic the- 
orem [3, [|| • Their generalization of FT is plausible, but 
their naming is misleading because both terms do not ac- 
curately reflect the entropy changes from non-adiabatic 
or adiabatic processes. 

AS na has an another name, the excess heat entropy 
change. But, as seen in [6|, the excess heat entropy 
change can be derived in adiabatic approximation. Prob- 
ably there is something wrong in 2] . In order to conform 
which one is an adiabatic contribution, one must per- 
form a rigorous verification in an adiabatic process (or in 
a non-adiabatic process), but they missed it. 

In this work, we study adiabatic and dibatic entropy 
changes along a single path. At first, the system is de- 
fined in Model System, and we will conjecture which one 
is the adiabatic path entropy change in Thermodynam- 
ics Interpretations. Then the conjectures are verified in 
Adiabatic Approximation by dividing the underlying pro- 
cesses into adiabatic and diabatic processes. In which, 
coarse graining and time scale separation [^, 0| will be 
appropriate mathematical tools for that. Lastly we con- 
clude with discussions about excess heat entropy change 
and adiabatic approximation in Conclusion. 

Model System.— Imagine a Markovian system of N 
possible states with an ergodic control parameter At. We 
assume that non-zero At gives the system work added 



(or done), W, and temperature of a reservoir (or reser- 
voirs) is constant. At any time t, the control parameter 
changes first from A t to At+ ro for unit time To, and then 
the system goes from previous state it to next state it+r 
by time evolution operator wj° + it (^t+T ) which is con- 
trolled by schedule of a control parameter. If the system 
undergoes a transition from a state to an another state, 
W it+T n (^*+ r o) pl a Y s a r °l e °f a transition probability. If 
the system is still in previous state, i.e. if it — it+r , then 
W i° +T0 i t (^t+T ) plays a role of a waiting probability. 

Let us consider a time evolution of the system which 
goes from an initial state io to a final state ix along a 
path 



'0 



(1) 



for a period of T = utq, where a is a pair of i and A. At 
the end of the schedule, the system goes back along the 
reversed path. For a single path, total entropy change is 
defined by 
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where is a function of the path [cr]^ @, Q. Here, 

P[(t]q and "P[ct]^ are respectively forward and reversed 
path probabilities, which are given by 
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where Pi o (0) and pi T (T) are respectively probabilities 
of the initial state io an d the final state ix- Since 
w it + i t (^t+r ) is not a transition probability but a time 
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evolution operator, there is no need to consider waiting 
time in path probabilities. 

Thermodynamics Interpretations.— Here, we will study 
adiabatic and diabatic path entropy changes in the view 
of thermodynamics. 

Total entropy change @ for a single path can be 
rewritten as 



AS tot [*}%= ^ In 



(5) 



As well known, it can be a sum of system and environ- 
ment entropy changes, AStot — AS sys +AS env [8|. Those 
are respectively given by 



AS sys [a]o = In 



and A^enuMo = / ] m 



PtT ( T ) 
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Reflecting on the first law of thermodynamics, AS = 
W + Q, the system can undergo a transition due to work 
W or heat Q. So, AS^^o-Jq can be regarded as a sum 
of two contributions, 



AS env [a]o = A5w[cr]o + AS Q [a 



(8) 



where ASVy[cr]J and ASqIctJq are respectively an adia- 
batic contribution caused by W and a diabatic contri- 
bution caused by Q (which are quitely different from Es- 
posito and Van den Broeck's adiabatic and non-adiabatic 
contributions) . 

According to the adiabatic theorem, the system be- 
comes an adiabatic process when At varies very slowly 
043 ■ But there is one more requirement, time scale sep- 
aration, t s . If the unit time or coarse graining period are 
larger than r s , the system goes to an adiabatic process 
[f|. If not, even At varies very slowly, the system is a dia- 
batic process, thereby there is a non-zero entropy change 
due to Q in non-equilibrium steady state. The non-zero 
entropy change is, as well known, the house-keeping heat 
entropy change [8l4l(| which is given by 



AS S [<= Vln 



m— 1 
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jtq (A mTo )pf mTQ (A mTo ) 



9) 



where pf(X) is steady state probability for A. In 
they have called it as an adiabatic contribution, AS a , but 
which is a contradiction in the view of thermodynamics, 
because ASq (or AS a ) is definitely due to Q. Addition- 
ally, if the system is coarse grained for a period r larger 
than t s , the system indeed goes to an adiabatic process 
as mentioned above, thereby ASq vanish jf|. Therefore, 
ASq (or AS a ) in non-equilibrium steady state can not 



be an adiabatic contribution at all, but rather a diabatic 
contribution. 

On the other hand, ASyv is definitely an adiabatic con- 
tribution, since non-zero At gives the system W as defined 
in Model System. From 0, © and ©, AS W is written 
as 



AS w [^]o 



(A, 



-Dm 



(A m ) 



(10) 



Here, ASyv is a type of excess heat entropy change ft^fllj]. 
but different from excess heat entropy change (that will 
be discussed in Conclusion). 

In summary up to here, total entropy change along a 
single path can be a sum of three different contributions, 

AStotMj = AS svs [o-]T + ASWMo + ASeMo- (11) 

If there is no house-keeping heat entropy change, the sys- 
tem is an adiabatic process, thereby path entropy change 
should be ASa — AS sys + ASyv- However, in a diabatic 
process, ASq is also to be considered. Therefore total 
path entropy change should be divided into an adiabatic 
path entropy change and a diabatic contribution, which 
are respectively 



AS A [o-]T =ln 
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(13) 

By comparing (|12[) and (fT3|) with [2j, it is clear that 
AS na = ASa and AS na = ASq. Therefore, AS na must 
be corrected to an adiabatic entropy change, and AS a is 
related to a diabatic contribution, if our conjectures, (|12[) 
and (Tl3j) . are right. 

Adiabatic Approximation.— Now, our conjectures 
above is verified by calculating path entropy changes in 
adiabatic approximation (or quasi-static approximation) . 

According to the adiabatic theorem, the system be- 
comes an adiabatic process when an external field varies 
very slowly It is well known that the adiabatic 

process corresponds to a quasi-static process in thermo- 
dynamics. Although the quasi-static process means a 
process to make the system staying extremely close to 
equilibrium steady states, it could be also applied to a 
non-equilibrium thermodynamics system So, per- 
forming calculations in the quasi-static process will give 
an adiabatic entropy change. 

At here, we have to deliberate on the adiabatic theorem 
in thermodynamics. As mentioned above in Thermody- 
namics Interpretations, in order to make the system an 
adiabatic process, the follows are to be required: Firstly 
At should be slowly changing, and secondly unit time or 
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coarse graining period have to be larger than the time 
scale separation. In fact, the second requirement con- 
tains the first one. So, if a process is coarse grained for 
a period r = uqTq larger than t s , the system becomes an 
adiabatic process @, 0] (we do not study in detail about 
t s in this work). Hence, when t > To, we can calculate 
an adiabatic entropy change exactly. The diabatic en- 
tropy change can be obtained as an additional part when 
t — To by comparing the total entropy change with the 
adiabatic entropy change. 

Let us introduce our strategy to calculate ASVv in 
a Markovian system. Firstly, introducing a temporar- 
ily time homogeneous Markov process, the time evolu- 
tion operator is coarse grained for a period r = tiqTq. 
Secondly, time evolution operators corresponding to re- 
spectively adiabatic and diabatic transitions, W s '(At) 
and N T (At), are defined. Lastly, an adiabatic entropy 
change, AS a, is calculated in adiabatic limit (or quasi- 
static limit). 

The first step starts from here. In the system, Markov 
chain at time t for unit time tq is given by 



P(t + T ) = W T °(A t+TD )P(t), 



(14) 



where P(t) = (pi(i), • • • ,p N {t)) Transpose is a column vec- 
tor of a state distribution, and W r °(A t ) = (w^(Xt)) is 
an N x N matrix of a time evolution operator. At this 
step, we consider a partial path 



H > H+t 



Am 



Ah 



> it+r , (15) 



for coarse graining period r = tiqTq along which (|T4|) is 
rewritten as 

P(t + r) = W T «(X t+T )---W T °(X t+To )P(t). (16) 

V 

no=T /t 

But it is rather complicated to calculate the matrices, 
because At varies in time. 

To make it simple, the temporarily time homogeneous 
Markov process is introduced. When the control param- 
eter is fixed at At+ r , the system will evolve along the 
following partial path 



lt+T_ 



H > H+t 



> it+r) ■ (17) 



In this partial path, (fTfij) is simplified as 

P(t + r) - W T (X t+T )P(t), (18) 

where 

W T (X t+T ) ee W^(X t+T ) ■ •■ W T °(A t+T ) = K(A t+T )) . 
v . ' 

no=T /tq 

(19) 

The coarse grained time evolution operator wJj(X t + T ) 
looks like a singular jump from a previous state to a 



next state, but obviously contains all transition and all 
waiting because it is a result of a path integration for all 
possible partial path, w^ (A t+T ) = J2[\] w Jj [A]t +T where 

WlM^ = W Z + T-T ' • • < +TO i( A *) ^ WhiCh tranSi " 

tions (or waitings) in the middle of the partial path (fl~7|) 
are just veiled. 

Here is the second step. If r — > oo or At = A = const, 
then the system heads to a unique steady sate and finally 
reaches to there regardless of an initial state P(0), 



P st (A) = W st (A)P(0). 



(20) 



where the steady state time evolution operator is 

W st (A) = W°°(A). (21) 

Here, P st (A) = (pf (A), • • • ,p%(\)) Tranapose is the 
unique steady state distribution for A. Since the equa- 
tion P0)| holds always, a trial distribution, P' = 
(• ■ • , 0,p'k = 1) 0, • • • ) T for any k, might be chosen in the 
place of P(0). Then, from $ZU$ and P', the following 
relation is derived, 



pf(X) = wf k (X) 



(22) 



for all i and k. For an any ergodic system, the relation 
is indeed valid and already known in mathematical fields 
of stochastic processes [12j • The relation means that all 
column vectors of W st (A) are same with P S '(A), 



W st (A) 



f pf(X) pf(A) 
Pi* (A t ) pi* (A) 



Pf (A) \ 

p!*(a) 



(23) 



V p%(X) p%(X) ■■■ p$(\)J 



To obtain this matrix, there is no need to multiply 
W T °(A) infinitely. That is simply obtained by the prin- 
cipal eigen-value equation, P st (A) = W r (A)P st (A). 

Now, let us consider a whole path from initial state io 
to final state it, 



0" 6 = *o 



At, 



Xt-t 



-t it 



(24) 



which is a sequence of the coarse grained partial paths 
(fl"7|) . where T = n'r. When t S> t or A t varies very 
slowly comparing to t, the system goes to adiabatic limit, 
and the time evolution operator becomes a steady state 
time evolution operator, 



W T (A t ) ~W st (A t 



(25) 



In this limit, the system is evolving near steady states 
of At, i.e quasi-static processes (in thermodynamics) or 
adiabatic process (in quantum mechanics). On the other 
hand, if r ~ To or At varies rapidly comparing to r, 
(|25p does not hold anymore, and we have to consider 
an additional part corresponding to a diabatic contribu- 
tion which is given by definition, N r (At) = W T (At) — 
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W s '(At). So, the time evolution matrix can be separated 
into two parts, 



W T (A t ) = W st (A f ) + N T (A t ). 



(26) 



This is the last step to calculate the path entropy 
changes. In adiabatic limit, i.e. when t > t or A t is 
very small comparing to t, total path entropy change is 
same with adiabatic path entropy change. So, from ([5]) 
and 1|22|) and ([25| . it is simply obtained as 



AS ad [a]T = In 



Pit( T ) 



E 



In 



Pi. 



( XjriT 



(27) 



On the other hand, in diabatic limit, i.e. if r ~ tq or At 
varies very rapidly comparing to r, we have to consider 
a diabatic contribution in entropy change which is given 
by definition, AS dl = AS tot - AS ad . From (|5l and (1271). 
the diabatic contribution is written as 



™ Wj j (ArnrM* n (A m ) 

m= 1 



By comparing (f2"T|). (f2"g)) with (IT2"1) . (fT5|) . we can ascertain 
that 



AS ad [< 



ASUHo and AS* [oft = AS Q [< (29) 



when t — tq. Therefore, our conjectures in Thermody- 
namics Interpretation, (fl~2j) and (fl~3j) . are right. 

Conclusion.— In this study, we conjectured which one 
is an adiabatic entropy change in Thermodynamics Inter- 
pretation, and then in Adiabatic Approximation we veri- 
fied our conjectures: AS a is an adiabatic entropy change, 
and ASq is a diabatic contribution. Therefore the terms 
of AS na and A5„ in Q should be completely switched 
as A5* na -> ASa and A5 a ->• ASq in fact. 

Even though the mathematical form of ASa is a type 
AS Seai d [Hi, ASU is definitely different from AS Qex 
because their origins are different from each other. In this 
work, because non-zero At gives the system only W, ASyv 
is not excess heat entropy change definitely but a type 
of excess heat entropy change. If non-zero At gives the 
system Q, ASq bx can be obtained as a diabatic entropy 
contribution, i.e. ASq — ASQ hk + ASq^. 

This work is consistent with other previous findings. 
One of them is M. Santillan and H. Qian's finding [6j. 
They obtained a similar result, S = — Q e x/T, in adia- 
batic approximation. In their study, they obtained the 
result in a molecular system with constant temperature, 
so that is a entropy change due to non-zero At and a type 
of excess heat entropy change. Moreover, T. Hatano's re- 
search in 1999 can be an evidence supporting our re- 
sults, too. He verified the Jarzynski equality L4J without 
consideration of house-keeping heat [l3(. Remembering 



that the Jarzynski equality (or the Crooks relation (ljj ) 
holds originally in adiabatic limit and not related with 
Q, T. Hatano's finding in [l| can be regarded as per- 
formed in adiabatic limit, and indeed performed in adia- 
batic limit (or quasi-static limit). Therefore, we conclude 
that ASa = AS sys + ASw is indeed an adiabatic entropy 
change, and ASq = A5g ox + ASQ hk is a diabatic con- 
tribution. 

Lastly, we faced with a interesting topic which is 
the time scale separation between equilibrium and non- 
equilibirum domains. Looking at (|25[) . we can ascertain 
that detailed balance condition is satisfied in adiabatic 
limit or in slow process which is already reported in [l7| . 
Researching about what determines r s can be an inter- 
esting topic. 
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